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Law of large numbers for branching symmetric Hunt 
processes with measure-valued branching rates 

Zhen-Qing Chen* * Yan-Xia Rent and Ting Yang 


Abstract 

We establish weak and strong law of large numbers for a class of branching symmetric 
Hunt processes with the branching rate being a smooth measure with respect to the underlying 
Hunt process, and the branching mechanism being general and state-dependent. Our work 
is motivated by recent work on strong law of large numbers for branching symmetric Markov 
processes by Chen-Shiozawa [6j and for branching diffusions by Englander-Harris-Kyprianou 
m- Our results can be applied to some interesting examples that are covered by neither of 
these papers. 
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1 Introduction 

1.1 Motivation 

The law of large numbers (LLN) has been the object of interest for measure-valued Markov processes 
including branching Markov processes and superprocesses. For branching Markov processes, the 
earliest work in this field dates back to 1970s when Watanabe [25, 26] studied the asymptotic 
properties of a branching symmetric diffusion, using a suitable Fourier analysis. Later, Asmussen 
and Hering [3J established an almost-sure limit theorem for a general supercritical branching process 
under some regularity conditions. Recently, there is a revived interest in this field using modern 
techniques such as Dirichlet form method, martingales and spine method. Chen and Shiozawa [6] 
(see also [23] ) used a Dirichlet form and spectral theory approach to obtain strong law of large 
numbers (SLLN) for branching symmetric Hunt processes. Among other assumptions, a spectral 
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gap condition was used to obtain a Poincare inequality that plays an important role in the proof of 
SLLN along lattice times. They proved SLLN holds for branching processes under the assumptions 
that the branching rate is given by a measure in Kato-class Koo(X) and the branching mechanism 
has bounded second moment. 

The spine method developed recently for measure-valued Markov processes is a powerful prob¬ 
abilistic tool in studying various properties of these processes; see, e.g., b mammalia usual- 
In mi, Englander, Harris and Kyprianou used spine method to obtain SLLN for branching (possi¬ 
bly non-symmetric) diffusions corresponding to the operator Lu + /3(u 2 — u ) on a domain D C 
(where j3 > 0 is non-trivial) under certain spectral conditions. They imposed a condition on how far 
particles can spread in space (see condition (iii) on page 282 of mi)- That the underlying process 
is a diffusion plays an important role in their argument and the branching rate there has to be a 
function rather than a measure. The approach of mi also involves p- th moment calculation with 
p > 1 which may not be valid for general branching mechanisms. Recently, Eckhoff, Kyprianou 
and Winkel [9] discussed the strong law of large numbers (SLLN) along lattice times for branching 
diffusions, which serves as the backbone or skeleton for super diffusions. It is proved in [9] Theo¬ 
rem 2.14] that, if the branching mechanism satisfies a p- th moment condition with p G (1,2], the 
underlying diffusion and the support of the branching diffusion satisfy conditions similar to that 
presented in mi, then SLLN along lattice times holds. 

In this paper, we combine the functional analytic methods used in [6] with spine techniques to 
study weak and strong laws of large numbers for branching symmetric Hunt processes as well as the 
L log L criteria. This approach allows us to obtain new results for a large class of branching Markov 
processes, for which (i) the underlying spatial motions are general symmetric Hunt processes, which 
can be discontinuous and may not be intrinsically ultracontractive; (ii) the branching rates are 
given by general smooth measures rather than functions or Kato class measures; (iii) the offspring 
distributions are only assumed to have bounded first moments with no assumption on their second 
moments. In addition, we use ^-approach instead of L p -approach for p G (1,2]. Now we describe 
the setting and main results of this paper in detail, followed by several examples illuminating the 
main results. 

1.2 Branching symmetric Hunt processes and assumptions 

Suppose we are given three initial ingredients: a Hunt process, a smooth measure and a branching 
mechanism. We introduce them one by one: 

• A Hunt process X: Suppose E is a locally compact separable metric space and Eg := LU{(1} 
is its one point compactification. m is a positive Radon measure on E with full support. 
Let X = (p.,'H,'Ht,6t,X t , n x ,£) be a m-symmetric Hunt process on E. Here {Ht : i > 0} 
is the minimal admissible filtration, {6t : t > 0} the time-shift operator of X satisfying 
XfO 6 S = Xt-\- s for s,t > 0, and £ := inf{t > 0 : Xt = d} the life time of X. Suppose for each 
t > 0, X has a symmetric transition density function p(t,x,y ) with respect to the measure 
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m. Let {Pt : t > 0} be the Markovian transition semigroup of X, i.e. 


P t f(x) :=U x [f(X t )\ = I f (y)p(t, x, y)m(dy) 

J E 

for any non-negative measurable function /. The symmetric Dirichlet form on L 2 (E,m ) 
generated by X will be denoted as (£,E): 


T = (u € L 2 (E,m) : lim - f (u(x) — Ptu(x)) u(x)m(dx) < +oo|, 
l t—>o tip J 


£(u, v) = lim - / (u(x) — Ptu(x)) v(x)m(dx), u,v€P. 
t ^° t Je 

It is known (cf. [5]) that (£,J~) is quasi-regular and hence is quasi-homeomorphic to a regular 
Dirichlet form on a locally compact separable metric space. In the sequel, we assume that X is 
m-irreducible in the sense that if A € B(E) has positive m- measure, then H x {Ta < +oo) > 0 
for all x € E, where Ta '■= inf{t > 0 : Xt € ^4} is the first hitting time of A. 


• A branching rate /i: Suppose [i is a positive smooth Radon measure on (E, B(E)). It uniquely 
determines a positive continuous additive functional (PCAF) A^ by the following Revuz 
formula: 


Here n m (-) 


Je 

f E U x (-)m(dx). 


lim —II 

t —>-0 t 


m 



f{X s )dA^ , 


/ € B + (E). 


• Offspring distributions {{p n (x) : n > 0},x € E}: Suppose {{p n (x) : n > 0},x € E} is a 
family of probability mass functions such that 0 < p n (x) < 1 and Yl'^LoPn( x ) = 1- F° r 
each x € E, { p n (x ) : n > 0} serves as the offspring distribution of a particle located at x. 
Let {A(x) : x € E} be a collection of random variables taking values in {0,1,2, •••} and 
distributed as P(H(x) = n) = p n (x). 


Define 

+oo 

Q(x) :=^2np n (x), xeE. (1.1) 

n =0 

Throughout this paper we assume that the offspring distribution { p n (x ) : n > 0} satisfies the 
following condition: 

Po(x) = 0, pi(x) 1 and supQ(x)<oo. (1.2) 

x£E 

From these ingredients we can build a branching Markov process according to the following 
recipe: under a probability measure P^, a particle starts from x € E and moves around in Eq like 
a copy of X. We use 0 to denote the original particle, X^{t) its position at time t and (j its fission 
time. We say that 0 splits at the rate p in the sense that 


P x (C@ > t\X<t(s) : s > 0) = exp (-A?). 
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When £0 > (, it dies at time (. On the other hand, when < (, it splits into a random 
number of children, the number being distributed as a copy of These children, starting 

from their point of creation, will move and reproduce independently in the same way as their 
parents. If a particle u is alive at time t, we refer to its location in E as X u (t). Therefore the 
tirne-t configuration is a E-valued point process X t = {X u (t) : u € Z t }, where Z t is the set 
of particles alive at time t. With abuse of notation, we can also regard Xf as a random point 
measure on E defined by X* := J2 u eZt Let (Ft)t >0 be the natural filtration of X and 

Foo = <?{X t : t > 0}. Hence it defines a branching symmetric Hunt process X = (17, J r 0O , Ft, Xt, P^) 
on E with the motion component X, the branching rate measure p and the branching mechanism 
function {p n (x) : n > 0}. When the branching rate measure p is absolutely continuous with respect 
to m , i.e., p(dy) = /3(y)m(dy) for some non-negative function /?, the corresponding PCAF A% is 
equal to (3(X s )ds , and given that a particle u is alive at t, its probability of splitting in (f, t + dt) 
is j3(X v (t))dt + o(dt). Since the function (5 determines the rate at which every particle splits, (5 is 
called the branching rate function in literature. 

Throughout this paper we use Bi{E) (respectively, B + (E)) to denote the space of bounded 
(respectively, non-negative) measurable functions on (E,B(E)). Any function f on E will be 
automatically extended to Eg by setting f(d) = 0. We use (f,g) to denote f E f(x)g(x)m(dx) 
and as a way of definition. For a, b € M, a A b := min{a,6}, a V b := rnaxja, b}, and 

log + a := log(a V 1). 

For every / € B + (E) and t > 0, define 


Mf) ■■= E /(*«(*))• 

uGZt 


We define the Feynman-Kac semigroup 1 ^ t by 


P t W-lV /(x) := n x leMAP^fiXt)] , f G B + (E). 


(1.3) 


Since X has a transition density function, it follows that for each t > 0, admits an integral 

kernel with respect to the measure m. We denote this kernel by p(‘3 _1 )^(t, x, y). The semigroup 
p(Q 1 asS ociates with a quadratic form (f^ -1 ^, J 7 ^), where P 11 = F n L 2 (E ; p) and 


£(<2 1 )^(u ) u) = £(u, u) — [ u(x) 2 (Q(x) — l)p(dx), u^F 11 . 

Je 

We say that a signed smooth measure v belongs to the Kato class K(X), if 

limsupn x [j4^] = 0. (1.4) 

40 xeE 

For every non-negative v € K(X), we have HGq^Hoo < 00 for every a > 0, where G a is the 
a-resolvent of X and G a v is the a-potential of v. Define £ a (u,u) := £(u,u) + a f E u 2 dm. By [23], 

/ u(x) 2 v(dx) < ||G a ^||oo^a(w, u), u^F. (1.5) 

Je 
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Thus when p is in K(A), F p = J and the quadratic form {£^ Oa 1 , jth'j } s bounded from below. 
For an arbitrary smooth measure p, we define 


Ai := inf ^ 1 - )M (ri, u) : u € F p with J u{x) 2 m(dx) = 1 


( 1 . 6 ) 


Assumption 1. Let p be a non-negative smooth measure on E so that the Schrodinger semigroup 
p(Q l)fi a j m j^ s a gy m metric kernel x,y) with respect to the measure m and is jointly 

continuous in (x,y) € E x E for every t > 0. Moreover, Ai € (—oo,0) and there is a positive 
continuous function h € F^ with J E h(x) 2 m(dx ) = 1 so that £ h) = Ai. 

Observe that if u is a minimizer for (USD , then so is |u|. Assumption 1 says that there is a 
minimizer for (11,61) that can be chosen to be positive and continuous. Clearly the following property 
holds for h: 

£(h,v) = I h(x)v(x)(Q(x) — l)p(dx) + Xi(h,v) for every v € F p . (1.7) 

Je 

The finiteness of Ai implies that the bilinear form F^) is bounded from below, and 

hence by j2j, { Pj^ : t > 0} is a strongly continuous semigroup on L 2 (E,m). Let a(£^ 

denote the spectrum of the self-adjoint operator associated with Let A 2 be the second 

bottom of a(5^ -1 ^), that is, 

A 2 := inf u) : uGF^, j u(x)h(x)m(dx) = 0, f u(x) 2 m(dx) = lj. 

J E J E 

Assumption 2. There is a positive spectral gap in A h '■= A 2 — Ai > 0. 

Define ^-transformed semigroup {P^F > 0} from { Pj® > 0} by 

Ait 


P t hf ( x ) = ^P t iQ ~ 1)fl (hf)(x) for x € E and / € B + {E). 


( 1 . 8 ) 


Then it is easy to see that {Pj 1 ', t > 0} is an m-symmetric semigroup, where fh := h?m , and 1 is 
an eigenfunction of Pj 1, with eigenvalue 1. Furthermore the spectrum of the infinitesimal generator 
of {PjjF > 0} in L 2 (E]fh) is the spectrum of the infinitesimal generator of { P 1 ^;t > 0} in 
L 2 (E]m ) shifted by Ai. Hence under Assumption 2, we have the following Poincare inequality: 


\\PtP\[L 2 (E,rh) < e Xht \\<P\\vHE,m) 
for all ip G L 2 (E, fh) with f E ip(x)fh(dx) = 0. 


(1.9) 


Remark 1.1. If the underlying process X satisfies that for each t > 0, the transition density 
function p(t, x, y ) is bounded and is continuous in x for every fixed y € E and that the branching 
rate measure p is in the Kato class K(X) of X, then it follows from [lj that the Feymann-Kac 
semigroup Pj'^ maps bounded functions to continuous functions and is bounded from L P (E , m) 
to L q (E,m) for any 1 < p < q < + 00 . By Friedrichs theorem, Assumptions 1 and 2 hold if we 
assume in addition that 


5 



the embedding of {F,£\) into L 2 (E,p) is compact. 

Such an assumption is imposed in [6] to ensure the spectral gap condition and to obtain Poincare 
inequality (11.91) . 


1.3 Main results 


Recall that o is the natural filtration of X. Observe that (cf. [23, Lemma 3.3]) for every 

x £ E and every / £ B + (E), 

E, [Mf)} = P?^f{x). (1.10) 

It is easy to see that Mf := e Alt Xf(/i) is a positive P x -martingale with respect to Ft- Let M aQ := 
lim t^+oo Mf. It is natural to ask when is non-degenerate, that is, when P x (M 0O > 0) > 0 for 
x € El Under the assumptions that (i) m(E) < oo, (ii) the Feymann-Kac semigroup Pj® 1 ' l>l is 
intrinsically ultracontractive, and (iii) h is bounded, it is proved in m that the condition 


+oo 


E JE 


Y kp k (y)h(y) 2 log + {kh(y))p(dy) 


< +oo 


k =0 


( 1 . 11 ) 


is necessary and sufficient for to be non-degenerate. Condition (11.111) is usually called the 
LlogL criteria. The first main result of this paper reveals that, in general, condition (11.121) below 
is sufficient for M 0 0 to be non-degenerate. 


Theorem 1.2. Suppose Assumptions 1-2 hold. If 

/ n +00 

h{y) 2 log + h(y)m(dy) + / Y k Pk{y)Kv) 2 ^og + (kh(y))p(dy) < +oo, (1.12) 

. J -Ie,. o 

then Mf converges to in L 1 (P X ) for every x £ E, and, consequently, P x (M oa > 0) > 0. 

Thus under condition (11.121) . Xf(/i) grows exponentially with rate — Ai. Note that when h is 
bounded, (I1.12P is equivalent to (11.111) . The next question to ask is that, for a general test function 
/ £ B + (E), what is the limiting behavior of X f (/) as t —>• oo? By (11.81) and (11.101) . it is not hard 
to deduce (see (I2.8j) below) that for every / £ B + (E) with f < ch for some constant c > 0, 

e Alt Ej, [Xf(/)] = h(x)P t h (f/h)(x ) h(x)(f , h) as t +oo. 


So, the mean of X t (/) also grows exponentially with rate — Ai. Our previous question is related to 
the question: for / £ B + (E) with f < ch for some constant c > 0, does Xf(/) grow exponentially 
with the same rate? If so, can one identify its limit? We first answer these questions in Theorem 
PI and Corollary 11.41 in terms of convergence in L 1 (P X ) and in probability, under an additional 
condition ([1.131) . 

Note that under Assumption 1, for every t > 0, Pf 1 has a symmetric continuous transition density 
function p h {t,x,y ) on E x E with respect to the measure m, which is related to p^~ l ^{t,x,y) by 
the following formula: 


P h (t,x,y ) 


c Ai t P iQ 1]ti (t,x,y) 
h(x)h(y) 


x,y € E. 
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Theorem 1.3 (Weak law of large numbers). Suppose Assumptions 1-2 and (11.121) hold. If there 
exists some to > 0 such that 


/ p^ Q l ^(to,y, y)m(dy) < +oo, or equivalently, / p h (t 0 ,y,y)m(dy)<+ oo, 
Je Je 

then for any x € E and any f € B + {E) with f < ch for some c> 0, we have 

lim e Alt X t (/) = M 00 (f, h) in L\F X ). 

>■+oo 

Corollary 1.4. Under the assumptions of Theorem \ 1.31 it holds that 

lim = probability with respect to P^, 

t->+oo E x \K t {f)\ h(x) 

for every x € E and every f £ B + (E) with f < ch for some c > 0. 


(1.13) 


For almost sure convergence result, we need a stronger condition (|1.14p below. 


Theorem 1.5 (Strong law of large numbers). Suppose Assumptions 1-2 and (11.121) hold. If there 
exists ti > 0 such that 


sup 

y&E 


p ( Q l )^{t^y iy ) 

h{y) 2 


< +oo, 


or, equivalently, supp h (t\, y, y) < +oo, 
y eE 


(1.14) 


then there exists flo C fl of P x -full probability for every x € E, such that, for every u € flo and 
every f € Bb(E ) with compact support whose set of discontinuous points has zero m-measure, we 
have 

lim e Alt X 4 (/)(w) = M 00 (u)(f,h). (1.15) 

»+oo 


Corollary 1.6. Suppose the assumptions of Theorem 11.51 hold and let fto be defined in Theorem 
11.51 Then 

lim X *(/)( u; ) = M °qM 

E x [X t (/)] h(x) 

for every oj G flo and for every f £ Bb(E) with compact support whose set of discontinuous points 
has zero m-measure. 


Remark 1.7. The condition (11.131) is equivalent to 

p h (to/2,x,y) 2 fh(dy)m(dx) < +oo. 


E JE 


Hence by [3 Page 156], P t h is a compact operator on L 2 (E,m) for every t > to/2. 
Assumption 2 is automatically satisfied if either (11.131) or (jl. 141) holds. 

To understand condition ([1.14D . we give some equivalent statements of (11.141) 
sumptions 1-2. 


Consequently 


under our As- 
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Proposition 1.8. Suppose Assumptions 1-2 hold. The following are equivalent to (11.141) . 

(i) There exists 1 . 1 > 0 such that for any t>t\, 

sup \p h (t, x, y) — 1| < c\e~ C2t (1.16) 

x,y^E 


for some C\,C 2 > 0. 


p h (t, x, y) —» 1, as 1 1 +oo uniformly in (x, y) € E x E. 


(1.17) 


(iii) There exist constants t,ct > 0 such that 

p(^ _1 ) M (t, x, y) < cth(x)h(y) for every x,y € E. 


(1.18) 


Property (11.181) is called asymptotically intrinsically ultracontractive (AIU) by Kaleta and 
Lorinczi in m- If the inequality in (I1.18P is true for every t > 0, and every x,y € E, then 
{P f (Q - l)l1 : t > 0} is called intrinsically ultracontractive (IU). It is shown in [15] that in case of 
symmetric a-stable processes (a € (0,2)), (AIU) is a weaker property than (IU). 


1.4 Examples 


In this subsection, we illustrate our main results by several concrete examples. For simplicity, we 
consider binary branching only, i.e., every particle gives birth to precisely two children, in which 
case Q(x) = 2 on E. Since o ^ Pk( x ) log + k = 2 log 2 on E, condition (jl.121) is reduced to 


L 


h{y) 2 log+ h(y) ( m(dy) + n{dy)) < +oo. 


(1.19) 


Example 1. [WLLN for branching OU processes with a quadratic branching rate 
function] Let E = M rf . In Example 10 of [11] . (A, II^ ) is an Ornstein-Ulenbeck (OU) process on 
M d with infinitesimal generator 

C = -a 2 A — cx ■ V on M d , 

2 

where a, c > 0. Without loss of generality, we assume a = 1. Let m(dx ) = (f )^ 2 e~ c ^ 2 dx. Then 
X is symmetric with respect to the probability measure m, and the Dirichlet form (£, T) of X on 
L 2 (M. d ,m ) is given by 

T = < / € L 2 (R d , m) : [ |V/(x)| 2 m(dx) < +oo 1 , 


£(u, u ) = 


1 


|V/(x)| 2 m(dx). 










Let /3(x) = b\x\ 2 + a with a, b > 0 be the branching rate function. Let P/ 3 be the corresponding 
Feynman-Kac semigroup, 

Ptf(x) ■= n x 

Suppose c > y/2b and a = \Jc 2 — 2b. Let 

A c := inf{A 6 M : there exists u > 0 such that (£ + /3 — A )u = 0 in R^} 


exp 


P(X 8 )ds f{X t 


be the generalized principal eigenvalue. Let 4> denote the corresponding ground state, i.e., 4> > 0 
such that (£. + /3 — X c )(j) = 0. As indicated in [11], A c = ^(c—a) + a > 0 and i f>(x) = (^) d ^ exp(A(c — 
a)|x| 2 ). Note that 4> € and (/> = e~' Kct Pf cf) on R rf . It is easy to see that in this example Ai = —A c 
and h(x) = 4>(x). The transformed process (Af^jIlJ) is also an OU process with infinitesimal 
generator 

L = -A — ax ■ V on M d . 

2 

Note that its invariant probability measure is fh(dx ) = h{x) 2 m{dx) = (—) d / 2 e~ a ^ x \ 2 dx. Let 
p h (t,x,y ) be the transition density of X h with respect to fh. It is known that 


p h {t,x,y) 



In particular, 

= ( ^g-2 at ) eX P (^F^l^l 2 ) • 

Thus f Rd p h (t, x, x)rh(dx) < +oo for t > 0. Moreover, we observe that condition (11.191) is satisfied 
for this example. Therefore Theorem 11.31 holds for this example. 

This example does not satisfy the assumptions in |6j. To be more specific, here the ground state 
h is unbounded and f3(x ) = 6|x| 2 + a is not in the Kato class Koo(A) of X. 


Example 2. [WLLN for branching Hunt processes with a bounded branching rate func¬ 
tion] Let E be a locally compact separable metric space and m a positive Radon measure on E with 
full support. Suppose the branching rate function /3 is a non-negative bounded function on E. Sup¬ 
pose the underlying Hunt process (X, n x ) satisfies that for every t > 0, there exists a family of jointly 
continuous, symmetric and positive kernels p(t, x, y) such that Ptf(x ) = f E p(t, x, y)f(y)m(dy), and 
that there exists si > 0 so that 



x, x)m{dx) < Too. 


In this case the Feyman-Kac semigroup 


Pff(x) := U x 


exp 


P(X s )ds /(A, 


( 1 . 20 ) 
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has a jointly continuous and positive kernel pP(t,x,y). It is easy to see that 

e -||0||oo t p(t, x , y) < pP(t, x, y ) < x, y ) for every t > 0 and x,y G E. (1.21) 

Properties fll.20[) and (11.211) imply that f E pP(s±,x, x)m{dx ) < +oo. Thus Pf is a compact operator 
on L 2 (E,m) for every t > s\. By Jentzch’s theorem (see, for example, | 22[ Theorem V.6.6]), — Ai 
is a simple eigenvalue of C + j3 where C is the infinitesimal operator of X, and an eigenfunction h 
of C + /3 associated with —Ai can be chosen to be positive and continuous on E. Suppose Ai < 0. 
We assume in addition that there exists S 2 > 0 such that 

/ p(s 2 , x, x) 2 m(dx) < +oo. (1-22) 

Je 

It follows from (|1.21l) and Holder’s inequality that for every t > S 2 , Pf is a bounded operator from 
L 2 (E,m ) to L 4 (-E,m). Thus h = e Xlt Pfh € L 4 (E,m) and so condition (|1.19l) is satisfied. Hence 
Theorem 11.31 holds. 

Conditions (11.201) and (jl.221) are satisfied by a large class of Hunt processes, which contains 
subordinated OU processes as special cases. By “subordinated OU process” we mean the process 
Xt = Ys t , where Yj is an OU process on M d and St is a subordinator on R + independent of Yj. 
In the special case St = t, X t reduces to the OU process. In general, the sample path of X t is 
discontinuous. Suppose the infinitesimal generator of Y t is 

jC = ^(t 2 A — bx ■ V on R d 

where a, b > 0 are constants, and St is a subordinator with positive drift coefficient a > 0. 

As is indicated in Example 1, Yt is symmetric with respect to the reference measure m(dx) := 
d /2 

(-T-) exp(— b\x\ 2 /(j 2 )dx. We use p(t,x,y) to denote the transition density of Y t with respect to 

m. It is known that 

P(t, x, y) = (l-e~ 2bt ^J exp ^ _ (ja;| 2 + \y\ 2 - 2x • ye bt ^ . 

By definition, the transition density of Xt with respect to m is given by 

P(t,x,y) =^\p(S t ,x,y)} . 

It is proved in EDI Example 1.1] that (11.201) and (11.221) hold for a subordinated OU process. 
Therefore, Theorem 11.31 holds for branching subordinated OU processes with a bounded branching 
rate function. 

Example 3. [LLN for branching diffusions on bounded domains with branching rate 
given by a Kato class measure] Suppose d > 3, E C M. d is a bounded C 1,1 domain (that is, the 
boundary of E can be locally characterized by C 1,1 functions) and m is the Lebesgue measure on 
E. Let 

1 

C = - ^2 di( a ijdj) 
i,j =i 


10 



















with a,ij(x) € C' 1 (M rf ) for every i. j = 1, ■ ■ ■ ,d. Suppose the matrix ( dij(x )) is symmetric and 
uniformly elliptic. It is known that there exists a symmetric diffusion process Y on W l with 
generator C. Let X be the killed process of Y upon E, i.e., 


Yt, t < t e , 
d, t > t e , 


where te '■= inf{i > 0 : Yt 0 E} and d is a cemetery state. Then X has a transition density 
function p E (t,x,y) which is jointly continuous in (x,y) and positive for every t > 0. The following 
two-sided estimates of pE(t,x,y) is established in [21, Theorem 2.1], extending an earlier result of 


x and the boundary of E. There exist positive constants c*, i = l,--- ,4, such that for every 
(t,x,y) € (0,1] x E x E, 

cif E (t,x,y)t~ d/2 e~ C2lx ~ yl2/t < p E (t,x,y) < c 3 f E (t, x, y)t~ d/2 e~ c ^ x ~ y ^ 2 . 


1 A 


&E{y) 

Vi 


where 5 E (x) denotes the distance between 


Q. Zhang. Let f E (t,x,y ) := ( 


We say that a signed smooth Radon measure v on U. d belongs to the Kato class K d, a (« € (0, 2]) if 


lim sup 

r l** iEK ' 1 


'\x-y\<r 


W\(dy) 
x — y\ d ~ a 


= 0 . 


(1.23) 


In fact K d a is the Kato class of the rotationally symmetric a-stable processes on R d . We assume 
the branching rate measure p is a non-negative Radon measure in K^ 2 - For any / € B + (E), let 


Pff(x) := E x [exp (A?)f(X t )}. 


Then P( l has a transition density p y E {t,x,y) which is jointly continuous in (x,y) and positive for 
every t > 0. It is shown in |16l Theorem 4.4] that there exist positive constants Cj, i = 5, ■ ■ ■ , 8, 
such that for every (t, x, y) £(0,l]xEx E, 

C 5 f E (t , X, y)t-d/2 e -<*\*-v\ a /t < Xj y) < C 7 f E ^ y)t -d/2 e -c s \ X -y\Vt_ (1 . 24) 

The infinitesimal generator of Pj 1 is (£ + p)\ E with zero Dirichlet boundary condition. It follows 
from Jentzch’s theorem that — Ai is a simple eigenvalue of (C + p)\ E and that an eigenfunction h 
associated with — Ai can be chosen to be positive with \\h\\ L 2^ Et d x ) = 1- Immediately, h is continuous 
on E by the dominated convergence theorem. We assume Ai < 0. Recall that E is bounded. Using 
the equation h = e^ 1 P y h and the estimates in (|1 .24jl . we get that for every x € E, 


c 9 (l A 5 e (x)) < h(x) < ci 0 (l A S E {x)) 


(1.25) 


for some positive constants c 9 , cio- Let 

hn \ e Xlt p^{t,x,y) 

PEV^^y) ■=— 7 / \i / \ — for x,yeE. 

h{x)h{y) 

Immediately condition (11.191) holds by the boundedness of h and condition (|1.14l) holds by (|1.24l) 
and (jl .2511 . Therefore both Theorem 11.31 and Theorem 11.51 hold for this example. 
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Example 4. [LLN for branching killed a-stable processes with a bounded branching 
rate function] Suppose d > 1, E = M rf , m is the Lebesgue measure on R d and a £ (0, 2). Suppose 
Y is a symmetric a-stable process on R d , and c(x ) is a non-negative function in <* ( a function 
q is said to be in K^ a , if the measure v(dx) := q(x)dx is in K-d,a where K d, a is defined in (|1.23l) i. 
Let X be the subprocess of Y such that for all / £ Bb(R d ), 


Ptf(x) 


E x [f{X t )\ = E x 


exp 


c(Y s )ds 



It is known that the infinitesimal generator of X is C, = A“/ 2 — c(x), where A 0 / 2 := — (—A)“/ 2 is 
the generator of a symmetric a-stable process. Let the branching rate function /3 be a non-negative 
bounded function on R d . Let V(x) := c(x) — j3(x). Clearly, |U| £ K^ a . For any / £ B + (R d ), let 

ft 


PffW ■= E a 


exp 


P(X s )ds f(X t 


Note that for every t > 0, Pt is bounded from L 1 (M d ,(ix) to L°°(R d , dx), and satisfies the strong 
Feller property. It follows from Cl that for very t > 0, Pj 3 is bounded from L pi (M d , dx) to L P2 (M d , dx) 
for any 1 < p\ < p 2 < -foo. Thus under our Assumption 1, the ground state h is a positive bounded 
continuous function on R d . The semigroup {Pj 3 : t > 0} is the Feynman-Kac semigroup with 
infinitesimal generator C /J = A“/ 2 — V. Assume in addition that V(x) = c(x) — /3(x) satisfies the 
following conditions: liminf| a ,|_ 5 , +00 V + (x)/ log |x| > 0 and V~ has compact support. Then by [T5l 
Theorem 5.5], the semigroup P t 8 is (AIU). Hence both Theorem 11.31 and Theorem 11.51 are true for 
this example. It is known from m that in case of symmetric a-stable processes, (AIU) is weaker 
than (IU). For instance, V(x) = c(x) — /3(x) with c(x) = log |x|l{| x |>#| and /? has compact support 
in -B(0, R) for some R > 1 is such an example. 


The rest of the paper is organized as follows. In Section 2, we review some facts on Girsanov 
transform and /i-transforms in the context of symmetric Markov processes, and prove Proposition 
11.81 Spine construction of branching processes is recalled in Section 3. We then present proof for 
the LlogL criteria, Theorem 11.21 in Section 4. Weak law of large numbers, Theorem 11.31 is proved 
in Section 5, while Theorem 11.51 on the strong law of large numbers will be proved in Section 6. 
The lower case constants ci,C 2 ,--- , will denote the generic constants used in this article, whose 
exact values are not important, and can change from one appearance to another. 


2 Preliminary 


Recall h £ is the minimizer in Assumption 1. Since h £ P, by Fukushima’s decomposition, we 
have for q.e. x £ E, n^-a.s. 

h(X t ) - h(X 0 ) = M t h + N t h , t > 0, 


where M h is a martingale additive functional of X having finite energy and A ] h is a continuous 
additive functional of X having zero energy. It follows from (11.71) and (131 Theorem 5.4.2] that IV/ 1 
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is of bounded variation, and 

N t h = —Ai f h{X s )ds - [ h{X s )dA ( f- 1 ^, Vi > 0. 

Jo Jo 

Following Section 2] (see also [6j Section 2]), we define a local martingale on the random time 
interval [0, £ p ) by 

Mt := l h(kl C ^ [ °’ Cp) ’ (2 ' 1} 

where Cp is the predictable part of the life time £ of X. Then the solution R t of the stochastic 
differential equation 

Rt = 1 + [ R s ~dM s , t £ [0, £ p ), (2-2) 

Jo 

is a positive local martingale on [0, ( p ) and hence a supermartingale. As a result, the formula 

= RtdH x on Rt PI {t < C} for x € E 

uniquely determines a family of subprobability measures {11^ : x £ E} on (Q,R). We denote X 
under {n^; : x £ E} by X h , that is 


IT 


f(X t h ) =U x [R t f(X t ):t <(] 


for any t > 0 and / £ B + {E). It follows from [H Theorem 2.6] that the process X h is an irreducible 
recurrent m-symmetric right Markov process, where fh(dy) := h(y) 2 m{dy). Note that by (12.11) . 
(12.21) and Dolean-Dade’s formula, 


R t = exp(M t -^(M c ) t ) 


0 <s<t 


h(X s ) 
h(X s _) 


exp 1 — 


h(X s ) 
h(X s _) 


^ , t £ [0, C p ), 


(2.3) 


where M c is the continuous martingale part of M. Applying Ito’s formula to log h(Xt), we obtain 
that for q.e. x £ E, IL^-a-s. on [0, £), 


1 


log h(X t )—\og h(Xo) = M t --(M c ) t +J2 (tog 
By (12.31) and (12.41) , we get 


h(X s ) h(X s ) - h(X s _) 




h{X s -) 


on [0, C)- 


— Ait—A] 


( Q - i)m 


(2.4) 


Therefore for any / £ B + (E), 


IF 


f(X t 


e Xlt 

h(x[ 


-If, 




h(X t )f(X t ) 


1 ^ //O 1 \ 

__p Q-^ { hf)(x) = P t h f(x). 
h[x) 


(2.5) 


Let (£ h ,J- h ) be the symmetric Dirichlet form on L 2 (E,fh) generated by X h . Then (12.51) says that 
the transition semigroup of X h is exactly the semigroup { P t h : t > 0} obtained from P^ 1 ' ll ' 
through Doob’s /i-transform. Consequently, / £ J~ h if and only if fh, £ and 

£ h (f,f) = £ (Q - 1 Hfh,fh)- A, [ f (x) 2 h(x) 2 m(dx). 

JE 
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In other words, <& h : f i-)- fh is an isometry from ( £ h ,E h ) onto 1171 and from 

L 2 (E , fh) onto L 2 (E , m). Let a(£ h ) denote the spectrum of the positive definite self-adjoint operator 
associated with £ h . We know from Theorem 2.6] that the constant function 1 belongs to E h , 
and £ h (\, 1) = 0. Hence 0 G &{£ h ) is a simple eigenvalue and 1 is the corresponding eigenfunction. 
Therefore 

A i : = inf u) : u € T h with J u{x) 2 fh{dx) = 1 j> = 0. 


Let A(j be the second bottom of a (£ h ), i.e. 


A 2 := inf u) : u G T h with j u(x)m(dx) = 0 and j u(x) 2 m(dx) = 1 

In view of the isometry $> h , we have A 2 = A 2 — Ai. So Assumption 2 is equivalent to assuming 
A£>0. 

Define 

at(x) := p h (t,x,x) for t > 0 and x G E. (2.6) 

Note that by dLSD for any x G E and t,s > 0, 



Pt+s<P(x) | = \P^PMx)\ 


< / p h (s,x,y)\P t h ip(y)\m{dy) 

Je 

< i^j E V h {s,x,y) 2 m{dy)^ \\Pt<p\\L?(E,m) 

< a 2 s(x) 1/2 e- Xht \\ip\\ L 2 ( E}fh) . 


(2.7) 


We use (/, 9 )L 2 (ffi) to denote f E f(x)g(x)fh(dx). For every g G L 2 (E,m) we can decompose g as 
g(x) = (g , f)T?(m) + v{ x ) with p(x) satisfying f E ip(x)fh(dx) = 0. Then by (12.71) . for any t > s > 0 
and any x G E, 


p t9{x)~ (g,l)u*(m)\ = \ p t^{x)\<e Xh{t s) a 2s {x) 1/2 \\p\\ L 2 {E ^ ) 

< e- x ^- s ^ 2s {x) l ' 2 \\g\\ L 2 {E ^. ( 2 . 8 ) 


Proof of Proposition f~Ol We only need to prove that (11.141) implies (|1.161) . For any / G L 2 (E,fh), 
define c/ := J E f(x)fh(dx). Immediately, we have for any t > 0, J E (f — Cf)fh(dx) = 0, and 
/ — Cf G L 2 (E,m). By (11.91) . we have 

\\ P tf - c /llL 2 (E,m) < e ~ Xht \\f\\ L 2 (E,m)- (2-9) 

Let t\ > 0 be the constant in (11.141) . By the semigroup property, for any t>t\/2 and x G E, 

P h (t,x,y) = [ p h (t 1 /2,x,z)p h (t-t 1 /2,y,z)m(dz) 

Je 

= p t tl/2 Uy), ( 2 . 10 ) 
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where f x (') •'= P h (t\/2,x, •). Note that by Holder’s inequality, 

\p h (t,x,y) - 1| = | f p h (t/2,x,z)p h (t/2,z,y)m(dz) - 1| 

J E 

= I [ (.P h {t/2,x,z) - \)(p h (t/2,z,y) - l)m{dz)\ 
Je 


< 


(p h (t/2, x, z) — 1 ) 2 m(dz) 


1/2 


(p W 2 , 2, y) - 1) m(dz) (2.11) 


1/2 


Note that Cf x = f E p h (ti/2,x,y)m(dy) = 1 and f E f 2 (y)m(dy) = (x). By (12.1011 and (12.91) . for 

any t > t\, 


r \ 1/2 

l^(p h (t/2,x,z) - l) 2 m(dz)J = || Pfc- tl)/2 f x - c fx \\ L 2 {Ejfh ) 

< f x \\ L2{Efh) 


= e- Xh{t - tl)/2 d tl (x) 1/2 . (2.12) 

Similarly, for any t > t\, 

U(p h (t/2,z,y) - 1 ) 2 m(dz)V /2 < e -A,(t-ii)/ 2 // ii(;(/ )i/ 2 . ( 2. 13) 

Combining (12.111) . ()2.121) and (|2.13p . we have for any t > t±, 

\p h (t,x,y) - 1| < (2.14) 

Therefore (11.141) implies that for any t > t\, (j 1.16 1 ) holds. □ 


3 Spine construction 

To establish the L 1 convergence of the martingale Mf, we apply the “spine” and change of measure 
techniques presented in m for branching diffusions to our branching Hunt processes. The notation 
used here is closely related to that used in [13] . It is known that the family structure of the particles 
in a branching Markov process can be well expressed by marked Galton-Watson trees (see, for 
example, mm and the references therein). Let T denote the space of all marked G-W trees. For 
a fixed r € 7~, all particles in t are labeled according to the Ulam-Harris convention, for example, 
0231 or 231 is the first child of the third child of the second child of the initial ancestor 0. Besides, 
each particle u € r has a mark (X u , a u , A u ) where X u : [b u ,(u) —>• Ed is the spatial location of u 
during its life span (b u is its birth time and £ u its fission time), a u = Cu ~ b u is the length of its life 
span, and A u denotes the number of its offspring. We use u < v to mean that u is an ancestor of v. 

Since in this paper every particle is assumed to give birth to at least one child, for each tree r, 
we can choose a distinguished genealogical path of descent from the initial ancestor 0. Such a line 
is called a spine and denoted as £ = {£o = 0, , ^ 2 , ■ ■ ■ }, where £, is the label of the ith spine node. 

Define node^(^) := u if u € £ and is alive at time t. We shall use {X t : t > 0} and {n* : t > 0} 
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respectively to denote the spatial path and the counting process of fission times along the spine. 
Let T denote the space of G-W trees with a distinguished spine. We introduce some fundamental 
filtrations that encapsulate different knowledge: 


Ft • — & {( 77 , X u , <X U , Au) , Cm — X u (s ), S € [ 6 U , f]), f € [ 6 U , Cm)} > 


^oo := V -Pt := cr{-^; t > 0}; 
t>o 


Ft := o' {F t , nodet(C)} , F 00 = \J Ft', 

t> o 

Gt '■= °{X S : s < t}, Goo '■= V @f, 

t> o 

Gt ■= ° {Gt, (node s (C) : s < f), (A u : u -< node t (C))} , Goo ■= \J Gt- 

t> o 

We assume P x is the measure on (f, Foo) such that the filtered probability space (f, F^, (F t ) t > 0 , P* 
is the canonical model for the branching Hunt process X described in the introduction. We know 
from [14| that the measure P x on (F, Foo) can be extended to the probability measure P x on 
(T,Foo) such that the nth spine node is uniformly chosen from the children of the (n — l)th spine 
node. 

For every t > 0, as in m, we define 


77 (f) := exp ^Ai t + a[ Q 


h{X t ) 
h{x) ’ 


Z(t) := e Xlt 


Mh) 

h(x) 


and 7?(i) := e Xlt TT A, 

h(x) 

M-<node t (?) 


Then 77(f), Z(t ) and 77(f) are positive P x -martingales with respect to the filtrations Gt, Ft and Ft, 
respectively. Moreover, both 77(f) and Z(t) are projections of 77(f) onto their filtrations: 


Z(t) = F x [r](t)\Ft\ , 77 (f) = P x [ 77 (f) \Gt\ for f > 0 . 


We call 77(f) the single-particle martingale and Z{t) the branching-particle martingale. As in [ 13 ], 
we define a new probability measure Q x by setting 

dQ x \p t = 77(f)dP x |^ t for f > 0. (3-1) 

This implies that 

dQx\r t = z {t)dF x \T t for t > 0. (3.2) 

The influence of the measure change (13.11) lies in the following three aspects: Firstly, under Q x 
the motion of the spine is biased by the martingale 77 (f). Secondly, the branching events along 
the spine occur at an accelerated rate. Finally, the number of children of the spine nodes is size- 
biased distributed, that is, for every spine node v located at x, A v is distributed according to the 
probability mass function {pk(x) := kpk(x)/Q(x ) : k = 0 , 1 , • • • }, while other (non-spine) nodes, 
once born, remain unaffected. More specifically, under measure 

(i) The spine’s spatial process X moves in E as a copy of (X h , n x ); 
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(ii) The branching events along the spine occur at an accelerated rate Jl(dx) := Q(x)h(x) 2 g,(dx). 
This implies that given Goo, for every t > 0, the number of fission times n* is a Poisson 
random variable with parameter At, where At is a PCAF of X having Revuz measure p. To 
emphasize this dependence, we also write Af for At- 

(iii) At the fission time of node v in the spine, the single spine particle is replaced by A v children, 
with A v being distributed according to the size-biased distribution {pk{X^ V J) : k = 0,1, • • • }. 
Each child is selected to be the next spine node with equal probabilities. 

(iv) Each of the remaining A v — 1 children gives rise to independent subtrees, which are not part 
of the spine and evolve as independent processes determined by the measure P shifted to their 
point and time of creation. 

For more details on martingale changes of measures for branching Hunt processes, see |14| and H2j. 


4 LlogL criteria 


In this section, we prove Theorem 11.21 It follows from [5], Theorem 4.1.1] that if v is a smooth 
measure of X h , then for every g £ B + (E) and t > 0, 


n- 


/ 

Jo 


g(x h s )d^s 


n g{y)v(dy)ds 

y 


(4.1) 


Since t is arbitrary, the monotone class theorem implies that for any f(s,x ) 
l € £> + [0, Too) and g € B + {E), 


n- 


If 


f( s ,Xg)dA v s 


fjn 


s,y)v{dy)ds. 


l(s)g(x) with 


(4.2) 


Note that for every / £ R + ([0, Too) x E), there exists a sequence of functions f n (s, x) = l n (s)g n (x) 
with l n £ £? + [0,Too) and g n £ B + (E) such that f n (s,x ) converges increasingly to f(s,x ) for all s 
and x. Thus by the monotone convergence theorem, (14.21) holds for all / £ £> + ([0, Too) x E). 


Proof of Theorem \l.A Since M t is a positive martingale, it suffices to prove that E^Mqo = h(x) for 
all x £ E. By [8] Theorem 5.3.3] (or [19] Lemma 3.1]), it suffices to show that if (11.121) holds, then 


Qa, ( lirnsup M na < Too 1=1 for every a > 0 and x £ E. (4-3) 

\N 9n—>-+oo / 

Recall that Goo contains all the information about the spine. We have the following spine decom¬ 
position for M t , 





Qx 

M t | Goo 

= Qx 





e A " £ h(X u (t)) 
uez t 


e Xlt h(X t ) T (A u -l)e x ^E^ J>(t-C„) Xt _ Cu (/i) 

U<£n t 
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(4.4) 


= e^h(X t ) + E (A u -l)e x ^h(X Cu ). 

U^n t 

Let Qt := &{X S : s < t} and Qq be the trivial er-field. It follows from the second Borel-Cantelli 
lemma (see, for example, [8l Section 5.4]) and the Markov property that for any a > 0 and M > 1, 

/ +oo 


(e x '™h{X na ) > M i.o.) = Q fh [y j ^(e x ^h{X na ) > M 


\n =1 
/ +oo 


G(n- 1)(T I — +°° 


( E n * (n _ 1)lT > M ) = +oo ) . (4.5) 


\n=l 


Recall that m is the invariant probability measure of (X,II^). Thus by Fubini’s theorem and 
Markov property, we have 


m 


+oo 


.n =1 


(n — l)a 


E^ [e Xinff h(X a )>M 


+00 


E n l (e AlTUT /i(X„*) > M 

n=l 
+oo ,, 

E / '*-{e A i 7 iCT /i(y)>M} m ( c ^y) 

n=l jE 

E ~,Aog + /i(2/) — log M 

E r ' 


ra=l 


-Ai 


> n<7 


< (-Ai) 1 / log+ h{y)m(dy) < oo. 


Therefore by (14.5p we have ( e Xina h(X na ) > M i.o.) = 0. Consequently, 


[ limsup e AlTlfT /i(X no -) < +oo j = 1. 

\N3n—>+oo / 


It is easy to check that the function x i->- Q x ^lim sup pj3n _ > ._|_ 00 e XincT h(X nrT ) < +oo) is an invariant 
function for (A,IlJ~). Recall that X has a transition density function with respect to rh. By [5j 
Theorem A.2.17], 

Qa, ( limsup e Xina h(X na ) < +oo j = 1 for every x € E. 

\N3n—>-+oo J 

Suppose e € (0, — Ai). For simplicity we use Q and A* to denote respectively the fission time and 
offspring number of the All spine node. 


+oo 


+oo 


+oo 


E e l( ' l Aih(X ( ; J — E e + E e i)^{Aih{X c .)>e e ti} 

2—1 2 — 1 2—1 

= : I+ 11. 


(4.6) 


Recall that Q 00 contains all the information of spine’s motion {Xt : t > 0} in E. For any set 
B € £>[0,+oo) x £>(Z + ), define N(B) := jj{* > 1 : (Q,Ai) € B}. Then conditioned on 
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N is a Poisson random measure on [0,+oo) x Z + with intensity dAg ft(l s )4(%), where 

Jl(dx) = Q(x)h(x) 2 n(dx). We have 


'+oo 

12 1 {A i h(X (i )>e^i} 

,i=l 


= Ir 

* * i 


p+oo + 00 


{kh(X s )>e^} dA s 


k =0 


Note that by (14.21) and our assumption (I1.12D . 


‘ + 00 

^2 1 {A i h{X ( .)>e^i} 
.1 = 1 


= IE 


r*+oo "I -00 


y ^Pk(tjs)l{ kh (x„)>eesydAV; 


k =0 


r-\-oo /* + 00 

/ / ^P^hlog+ikHy^esjPidy)^ 

jE k =0 

.. +oo 

1 / '22 k Pk{y)^og + (kh(y))h{y) 2 n{dy) < +°o. 


= £ 


k =0 


Thus 


^ l {A ih {X Ci )>e<i} < + 00 ) = 1 fOT m - a - e ' XeE - 

This implies that II is the sum of a finite many terms and so Q X (II < Too) = 1 for fh- a.e. x £ E. 
On the other hand, since Q(x) is bounded on E, we have 

r+OO +°° 

e x ' s ^p k (X s )kh(X s )l {kh{ t s) < e es } dAt 
k= 0 
+oo 

k= 0 

p+oo r ~J~°o 

/ f \ l -Lff^O X ^ /V / \ v'-'V / \ 1 / \ 9 / 7 S ~jj§ 

k =0 


QmW = n 


h 

m 

r+oo 


/0 


< 


< 


/■+oo /■ ' y 

/ / eAl *X ]pk{y) k Ky) l {kh{ y )<e^}p{ d y) ds 

Jo Je k=0 

/»+oo /• ~l~°^ 

/ / e iXl+£)sS 22p k {y)Q{y)h{y) 2 p{dy)ds 

Jo Je fc=0 

IIQIloo [ h(y) 2 y(dy) 

Je 


< Too. 


Thus we have (Qh;(/ < Too) = 1 for m-a.e. x £ E. Now we have proved that 
Q x e Al ^ J 4j/i(X^.) < Too^ = 1 for m-a.e. x £ E. 

It is easy to check that the function x it- Aih{X^.) < Too^ is an invariant function. 

Thus by [5J Theorem A.2.17], we get Q x e Al ^ A^Xq) < Too^ = 1 for every x £ E. By 

(14.41) we have 


( 


lirnsup ( 

\N 3n—^H-oo 


M na Goo 


< Too ) = 1 for every x £ E. 


By Fatou’s lemma, we get Q x (lirri infpj 3 „^ +C o Af n(J < Too) = 1. Note that M n J is a positive 
Q x -martingale with respect to E na , so it converges almost surely as n T oo. It follows then 

( lirnsup M na = liminf M na < Too 1=1. 

\N9n—>+oo NOn->-+oo J 
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□ 


5 Weak law of large numbers 

In this section, we present a proof for Theorem 11.31 

Lemma 5.1. If Assumption 1 and (11.121) hold, then for every t > 0 and <j> € B^(E), 

/ E x [X t (,4>h) log + Xt(0/i)] h(x)m(dx) < +oo. 

J E 

Proof. First we note that for every x € E and (f> € B^f (E), 

E x [X t (</>h)log + X t (<f>h)] < e~ Xlt \\(/)\\ 00 h(x)¥ lx ^e Xlt h(x)~ 1 X t (h) log + X t ((j)h) 

= e~ Alt ||0||oo/i(z)Qx [log + Xi(>h)] . 

Recall that Goo contains all the information about the spine. Since for any a,b > 0 

log + (a + b) < log + a + log + b + log 2, 

we have by Jensen’s inequality 


[log+X t (#)] = 


< 


< Qx 


< 


— \£.X 


log + X t (4>h) | Q, 
log Qx X* (</>/&) V 11 
log + (Qx ^t((f>h) | Go, 
log + Qx X t ((j>h)\ Gc 


+ r 

+ log 2. 


By (15.21) and (|5.4|) . it suffices to show that 

/ Q x log + Qx X t ((f)h) I Goo h(x) 2 m(dx) < +oo. 


We get the spine decomposition for Qx X t ((j)h) I Gt 


as follows: 


M^h) | 0oo] = m{X t ) +^ A n- l)% Cu [X t - Cu (#)] 

u-<£t 

< m(x t ) + ]T M<i>\\ooe- Xl( t- M h(XcJ 

U<£t 


(5.1) 


(5.2) 


(5.3) 


(5.4) 


(5.5) 


(5.6) 


Note that log + (o6) < log + a + log + b for any a, b > 0. Using this and an analogy of (15.31) as well as 
the assumption that Ai < 0, we have 


log + Qx X t ((j>h)\Go 


< log + j <j>h(X t ) + ^2 


Aj (t—Ci 


MXc 


20 
































< log +(<f>h(X t )) + Y lo S + (^ull^llooe Al(i Cu) h(X (u )j + log + n t 

< log + {(t>h(X t )) + Y ( lo g + Il0l|oo - Ai (t - Cu) + log + (A u h(X Cu j)) + n t 

U<(,t 

< log + ((fh(X t )) + Y lo § + (AxM^cj) + (! + log + Halloo - Ai t)n t . 


(5.7) 




By (11.121) and using the fact that m(dy ) is an invariant distribution, 


L 


log + (4>h(Xt)) h(x) 2 m(dx ) = / log + (f>h)(x)m(dx) < + 00 . 

J Je 


Hence (15.51) is implied by 


L 


Y log + (H u /i(X c J) +n t 


h(x) 2 m(dx ) < + 00 . 


(5.8) 


Recall that conditioned on N(-) = j){i > 1 : (Q,Ai) € •} is a Poisson random measure on 
[0,+ 00 ) x Z+ with intensity dAf Ylk£Z + Pk(X s )5k{dy)- We have 


L 

L 


Y tog + (A u h(X c J) +n t 

u-<£t 


h(x) 2 m(dx) 


Y log + {A u h{^u)) +n t 


n 


h 


u^t 
rt +°° 


h(x) 2 m(dx) 


IE 

ft r +°° 


/ Yp^{ log + (kh(X s )) + l)dA>f 

"'° fc =0 


m(dx) 


< t 


[ ds f V'pfc(y) (log + (M( 2 /)) + l) Q(y)h(y) 2 Li(dy) 

Jo Je k=0 

' <• +00 ,, \ 

1 °g + (^(s/))^(2/) 2 A fc (‘^2/) + IIQ||°o j^h(y) 2 p{dy)j . 


□ 


Immediately the last term is finite by () 1.12 1) . Hence we complete the proof. 

Lemma 5.2. If Assumptions 1-2 and (11.121) hold, then for any s, a > 0, m G N, and any x € E, 
lim e XAs+t ^X s+t {(fh) - E x Ui( s +*)x s+i (#) I X t ] = 0 in L\F X ), (5.9) 


t—y 00 


lim e Ai(n+m)a x {(f)h) _ Ex r e Ai(n+m)<7 X (^) | ^ 

N 9 n— >-+oo v ' L \ y 1 

/or every f> € (E). 


= 0 P T -a.s. 


(5.10) 
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Proof. For any particle u € Z s , let {X“’ s : t > s} denote the branching Markov process initiated by 
u at time s. It is known that conditioned on J- s , X u,s and X v,s are independent for every u,v € Z s 
with u 7 ^ v. For every s,t > 0, define 


Ss,t ■■= e Al % +t (#0 = e Al * Y KitiM’ 

uGZt 


and 


S s ,t ■— e 1 Y, 


*}■ 


u&Zt 


Obviously S s j > S S} t- 

First by the conditional independence and the Markov property we have 


E, 


5 'mar,na [Sma ,na | Fn 


— IE iT Vox Smcr,n<j | na 


= e 2Xina E x Y Var ^(n+m)cr 


-U£Zn 


J-'n. 


< e 2Xina E 


= e 2 Ainff E 7 


X E 


.uGZ n 


•KjrliflLO 

x \ ( n-\-m)cr 




T 

J n 


Y, (^niff(^)) •*-{X mo .((/.h)<e- A 1" CT } 

-UOzZyuj 


Let g m cr,ncr(y ) IEy (^mcr( 4>h)) l{X mCT (0/i)<e _;, ' 1 no '} 


. Immediately we have 


(5.11) 


W,na(y) < e- Ainff E y [X mff (^)] < e - A iO+-)-|| 0 || oo/i ( 2/ ). 


Thus g m(7 

,na e L 2 (£ , m) , and 


11 9m<r,na \\L 2 (E,m) X e A iO+ m ) o '||0|| oo . 

Using (I2.8|) . we continue the estimates in (15.lip to get that for n € N with na > 1, 


Note that 

+OC 


< 


E, 


Sma,ncr E x 


5, 


ma.na J na 


T 

J n 


< e 2Xina E x [X nrT ( 9rn<J , na )} 

= e Xina h(x)P% a (g ma , na /h)(x) 

e Xina h{x) ((g m <T,na,h) + e Xh/2 a 1 (x) 1/2 e~ Xhncr \\g rnr7}na \\ L 2 ^ Etm) 


< e 


Ai na 


h(x)(g v 


, h) + e Xh ^ 2 h{x)ai{x) l ^ 2 e~ Xhna ~ Xl1 


+oo 


Y eXinrT (g™,n<r,h) = Y 1 


,Ai na 


E„ 


n =1 


n= 1 


(Xmcr(4>h)) 2 l{x ma (^)<e-Ai-} h{y)m{dy) 


(5.12) 


(5.13) 
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< 


e Ai(*-l)ff E 


j h(y)m(dy) ^ 
h{y)m{dy) 

r r+o o /• 

/ h(y)m(dy) / z 2 P. y (X m(7 (#) € dz) 

Je Jo J i 


+0O 2 


l{X m!T (</i/i)<e- A i SCT } 

(X m a(4>h)) l{X mCT (0/i)<ie- x l a } 

+oo ^ 


ds 

dx 


< 


1 


— Aicr 


< 


— Al<7 

1 

—Aict 

1 

— Ai<7 


' lVze A l CT — AiCJX 2 

r+oo 

h(y)m(dy) z¥ y (X ma ((j)h) £ dz) 

Jo 


dx 


—Aicr 


JE 

— X\(T—\\m(T 


h(y)E y [X ma +h)] m(dy) 


||0||oo / Hy) 2 rn(dy) < + 00 , 


(5.14) 


where in the second equality above we used the change of variables x = e 1 ^ <T . It follows from 

' + °° 117 (S m(7tnC j - E x [S matna \ JF n 


(15.1311 and (15.141) that X) n =L ^ 
lemma, 


< + 00 . Thus by the Borel-Cantelli 


lim S ma ncr IE;/; 

n—>-+oo 


J m<j,n<j J ncr 


T 

J n 


= 0 P T -a.s. 


(5.15) 


Note that for every n, m £ N, we have 


E., 


Smo,na 5 'mcr,na 


= E. r 


= e Xina E. 


E x[Sma ,na | JF na ] — E x [S ma ,ncr | ncr] 

x m(7 {4>h) l{x mCT (<t>h)>e- x i na } 


y ] ^X u (ncr) 

-U^zEficr 


(5.16) 


Let fma,no{y) E y ~Km<r(4 > h)J^Xma(<l>h) > e- x i n!T } 


. Obviously, 


frnrr/na {y) < ll<^llooE y [x mcr (h)} < e 


— Aimer | 


Mv) 


(5.17) 


and so f m(T ,na £ L 2 (E,m). Then by (12.811 . we have for n £ N with na > 1, 

RHS Of = e X ^E x [X nc7 (fma,na)] = h{x)P+(f ma , na /h){ X ) 

< h(x){fma,na,h) + e Xh,2 d l {x) l,2 h(x)e~ Xhna \\f rna)na \\ L 2 {E ^ m) 

< h( x )(fma,na,h) + e Xh l 2 dx (x) 1/2 h(x)e~ Xhna ~ Xima || (j)\\ oo, 

where in the last inequality we used (|5.17l) . Note that 


(5.18) 


+00 


+oo 


^ fma, not h) — ^ ^ ^ma {(frh) l{X mcr ((f)h)>e~ x i ncT } h(y) rn (dy) 


72—1 


72—1 ' 


r+oo 


< 


ds / E ? , 


X TOCT (^/i)l {Xro<rW>e -A 1CTS} J h(y)m(dy) 

-|-oo r+oo 


r r -|-oo r+oo 

/ h(y)m(dy) / ds z¥ y (X ma ((j)h) £ dz) 

Je Jo Je~ x i cts 

r r+oo ^ r+oo 

/ h(y)m(dy) / —-—-dt / z¥ y (X mrT ((j)h) £ dz) 

J E J 1 1 J t 
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Y r r+oo rz ^ 

J h(y)m(dy) J z¥ y (X m<T (0/i) € dz) J -dt 

^ r p+oo 


—X\a 

1 


—^1°" Je 

The last term is finite by Lemma 15.11 Thus we get 

+oo +oo 


r r- too 

J h(y)m(dy) J z log + z¥ y (X ma ((j)h) € dz) 

E y [X mr7 (c/)h) log + 

^ma (#-)] h(y)m(dy). 


^ ®i[* ^ma,na ^rna.nrr} — ^ ] E x E x [S mrr _ n(T | F no-] E x[Smcr,na \ F no-] 


< +OO. 


(5.19) 


n =1 


n =1 


Recall that S mr7tna > S m(Tjna for every n > 0. Again by the Borel-Cantelli lemma, we have 

lim S ma ncr S mcr ncr — lim E x\Sma,na I F ncr ] ¥j x \S ma nrJ I F ncr ] — 0 P x -a.s. (5.20) 

n—>+oo n—H-oo ’ 1 ’ 1 


Now (|5.10p follows from (15.15(1 and (15.20j) . 

Substituting na by t and ma by s in (|5.13l) . we get for any t > 1, 


E, 


S s ,t ~ E* 


S s ,t Ft 


< e Xlt h(x){g s , t ,h) + e Xh/2 h{x)Mx) 1/2 e~ Xht ~ Xls ' 


Using a similar calculation as in ([5.141) . we get for s > 0, f 1 +<x> e Xlt (g s j, h)dt < Too. Consequently 
limj^+oo e Xlt (g s , t , h) = 0. Hence E x (s S:t - E x S s , t \ F t 


0 as i —» Too, or equivalently 


lim S 8it - E x [S Stt \F t ] = 0 in L 2 (P x ). 

t —^-|-oo V J 

Substituting na by t and ma by s in (j5.18|) . we get for any t > 1, 


(5.21) 


E„ 


S s ,t - S a>t 


= E, 


E x [S s ,t | Ft] - E x 


S a ,t | F t 


< h(x)(f St t, h) + e Xh ^ 2 ai(x) 1/2 h(x)e~ Xht ~ XlS ' 


By similar calculation as in (15.191) . we get / 0 + °°(/.s,t> h)dt < Too for all s > 0. Hence lim t-^+oo(fs,t, h) 
0. Thus we have 


E, 


S s ,t - S a>t 


= E, 


E x [S s ,t | F t ] — E x 


S s ,t | Ft 


0, as t —>• Too. 


Therefore (15.91) follows from (15.211) and (j5.22[) . 

Lemma 5.3. Under the conditions of Theorem 1 1 . 31 we have 


(5.22) 

□ 


lim lim e Xl ^K x \X s+t ((fh) I F t ] = M^h, h) in L l { P*), 

s—>-+oo £—>-+oo 


(5.23) 


for every cf € (E) and every x £ E. 
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Proof. Recall that M £ converges to M 00 in L 1 (P X ) by Theorem 11.21 It suffices to prove that 


lim lim E^ 

s— S-+ 0 O t —>+oo 

Note that by the Markov property, 


e Al(s+t) E x \K a+t (</>h) I F t ] - M t ((j>h, h) 


= 0. 


(5.24) 


e^+^E* [X s+t (#) \F t ] = e x ^ s+t ^E XAt) [X s ^h)} 

uGZt 


= e Al *X* (hPg 


(5.25) 


Thus we have for any s,t > to, 


E a 

< e Alt E, 

< e Al4 E 


e Al(s+4) E x [X s+t ((fh) | F t ] - M t (cj)h, h) 
t (hPjff) — h((j)h, h 
hPg(j) - h{4>h, h) 


X £ 


Using 


E„ 


= h(x)Pt{\p^-(<l>h,h)\)(x). 
we continue the estimation above to get: 
e Al(s+4) E x [X s+t (cfh) | Ft] - M t ((f>h, h ) 


< 


e -^.-t«rn m(xh{x) pH {z ^ ){x) 


< e ~ x h(s~to/2)\ 


M x ) 


at 0 (y) 1/2 m{dy) + a to (x) 1/2 e Xh(t to/2) ^a to (y)m(dy)J 


\ 1/2' 


0 as t —> Too, and then s — >• Too. 


Proof of Theorem Let (f> = f /h, then <f> € Bt ( E ). Note that for any s, t > 0 


□ 


e^^Xt+styh) - 


< 


T 


e A i( t+s )l t+s (^) - e A i( s + t )]E a . [X s+t (#) | J) 
e Al(s+t) E x [X 8+t (<f>h) | Ft] ~ M^h, h) 


Therefore Theorem 11.31 follows immediately from Lemma 15.21 and Lemma 15.31 


□ 


6 Strong law of large numbers 

In this section, we prove Theorem 11,51 
6.1 SLLN along lattice times 

Lemma 6.1. Suppose the assumptions of Theorem \1.5\ hold. Then for any a > 0 and any x € E, 

lim e Xina X na ((j)h) = Moo((j)h, h) ¥ x -a.s. (6-1) 

n—>-+oo 

for every cj) € B^ ( E) 
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Proof. Recall that M, x = lim n ^ > . +00 e Xina X na (h). Let g(y) := (f)h)(y) — h(y)((j)h,h). The conver¬ 
gence in (16.111 is equivalent to 


lim e Xina X na (g) = 0 P x -a.s. 

n—>-+oo 

For an arbitrary m € N, 

g Ai ( n + m )°'X (n+m)o . (g) 

= (e Al(n+m)CT X (n+m)o . (g) - E x [e x ^ n+m ^X {n+m)a {g) | T no \ ) + E, 

= : In T I In- 


( 6 . 2 ) 


e x ^+™)°X {n+m)a (g) 


F n 

(6.3) 


Note that by Lemma 15.21 we have 
I n = e x ^ n+m ^X {n+m)a ^h) -E x [e X ^ n+m) ^ {n+m) Mh) \ F n 
On the other hand, by Markov property, we have 

Iln = E e Al(n+m)CT E^ (nff) [X ma ( 9 )] 

UlziZncr 

= E e Al(n+m)(T P^- 1)M 5(^(^)) 

U^Zficr 

= e x ^ n+m >X nc {P$-^g). 


0 as n t Too P x -a.s. (6.4) 


(6.5) 


Note that our assumption (11.141) implies ()1. 17f) . Then for any fixed e > 0, there exist m > 0 
sufficiently large such that 

sup \p h (mcr, x, y) — 1| < e. 

x,y£E 

Then for any x € E, 

\ e Mm* p{ Q-l)» g{x) \ = Ih^ph^g/h)^ 

= h(x) 


(p h (ma,x,y) - l) 4>(y)m(dy) 

< h(x) p h (ma,x,y) - 1 ct>(y)h(y) 2 m(dy) 

Je 

< eh(x)((j)h,h). 


Consequently by (16.511 . we have 


\II n | < e{<j>h, h)M na . 

It follows from (16.31) . (16.41) and (16.61) that 

lim sup |e Aincr X ncr (g)| < £(4>h,ti)M 00 P x -a.s. 

n—H-oo 


( 6 . 6 ) 


Hence we get (16.21) by letting e —> 0. 


□ 
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6.2 Transition from lattice times to continuous time 

In this subsection we extend the convergence along lattice times in Lemma fG.ll to convergence along 
continuous time and give a sketch of proof of Theorem II.51 The main approach in this subsection is 
similar to that of |6 S Theorem 3.7] (see also [3j Theorem 1’]). According to the proof of (6j Theorem 
3.7], to prove Theorem 11.51 it suffices to prove the following lemma: 

Lemma 6.2. Under the conditions of Theorem, \ 1.51 for every open subset U in E and every x € E, 

liminf e Ali X t (l^//i) > Moo [ h,(y) 2 m(dy) ¥ x -a.s., (6-7) 

t~>+ o° Ju 

Note that if (16.71) is true for any bounded open set U in E with U C E, then (16.71) is true for 
any open set U. In fact, for an arbitrary open set U in E, there exists a sequence of bounded open 
sets {U n : n > 0} such that U n C E and U n t U. Hence if U n satisfies (16.71) . we can deduce that U 
satisfies (16.71) by monotone convergence theorem. 

In the following we assume that U is an arbitrary bounded open set in E with U C E. For any 
e, a > 0, n € N and x € E, define 

U £ (x) :={y eU : h(y) > K x )}i 

^ n,u := i _|_ £ (hlu)(X u (na))l{ Xv (t)eU e (x u (na)) for all vecf ,na and t£[na,(n+i)a]}i 

and 

ocr,£ . \ A Ai na-ya,e 

°n ‘ / y e * 71 , 11 * 

Here for each t > na, £“ ,n<T denotes the set of particles which are alive at t and are descendants of 
the particle u € C na . 

To prove Lemma 16.21 we need the following lemma. 

Lemma 6.3. Suppose that U is an arbitrary bounded open set U in E with U C E. Under the 
conditions of Theorem 1 1. 51 we have 


lim S°’ e - E x [S^ e I T na \ = 0 P x -a.s. (6.8) 

n->+oo L 1 J 

for every x € E. 


Proof. Note that for n € N and {Yj : i = 1, ■■■ ,n} independent real-valued centered random 
variables, 

E| = Var = ]Te|1-| 2 . 

i— 1 i —1 i =1 

Thus by the conditional independence between subtrees and the Markov property of a branching 
Markov process, we have 




{SZ £ - Ex [5^1^])' 


En 


g 2Xincr 


E 


e t 


\Y a ’ e — \V a ’ £ 

I L n,u [ 2 n,u 


T 

J n 


E n 
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< 2e 2Xinrj 

E 

E \\Y a,£ \ 2 1 T 1 

[1 1 71,U 1 | ** n<r\ 




= 2 e 2Xina 

E 






< 2(1 + e 

r 2 e 2Mna y ^^(X^na)) 



7L(zZ n cr 


Thus 


+OO 


X>x {S^~E x [S^\T na ]y 


n= 1 


22 ( h 2 lu){X u {na )) 

JlL^Zncr 


< 2(1 +£)- 2 ^ e Al2n<T E : 

n=1 
oo 

< 2(1 + e)- 2 ^ e Xl2na P^- l ^{h 2 lu){x) 

n—1 

oo 

< 2(1 + e)~ 2 h(x) ^2 e X ' 2n °PUhlu)(x). 

n= 1 


By Proposition 11.81 for any e > 0, \p h (ncr, x, y) — 1| < e for every x,y € E when n is sufficiently 
large. It follows that P^ a (hlu)(x) < (1 + e) f u h(y)fh(dy) < +oo, and so 


+oo 


(S°’ e — E x [S% £ | Pna]) <ch( X ) 


OO 

22 e XincT < oo, 


71= 1 71=1 

where c is a positive constant. This together with the Borel-Cantelli lemma yields (16.8p . 


□ 


Proof of Lemma 16.31 : If U is an arbitrary bounded open set in E with U C E, then (16.71) follows 
from Lemma 16.31 in the same way as [6] Theorem3.9]. We omit the details here. By the argument 
after UEZD, we know that (16.71) holds for any open subset U of E. 
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